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Abstract 

In this paper we study the semileptonic decays of B~ — > (tj c , J/ip)l~vi. We firstly evaluate 
the B c — > (rj c ,J/^) transition form factors Fo(q 2 ), F + (q 2 ), V(q 2 ) and ^4o,i,2(<? 2 ) by employing 
the pQCD factorization approach, and then we calculate the branching ratios for all considered 
semileptonic decays. Based on the numerical results and the phenomenological analysis, we find 
that: (a) the pQCD predictions for the values of the B c — > r\ c and B c — > J/^f transition form 
factors agree well with those obtained by using other methods; (b) the pQCD predictions for 
the branching ratios of the considered decays are Br (B~ — > i] c e~ i> e (u~ v^)) = (4.41^'qq) x 1CP 3 , 
Br{B- ^n c T~u T ) = (1.37±g;||) x 1(T 3 , Br(B- J/* e -p e (/i-^)) = (10.03±^;f|) x 1(T 3 , and 
Br(B~ J/^t~v t ) = (2.92+°;^) x HT 3 ; and (c) we also define and calculate two ratios of 
the branching ratios R r)c and Rj/^,, which will be tested by LHCb and the forthcoming Super-B 
experiments. 
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I. INTRODUCTION 



The charmed B c meson was found by the CDF Collaboration at Tevatron |l( fifteen 
years ago. It is a pseudoscalar ground state of two heavy quarks b and c, which can decay 
individually. Being below the B — D threshold, B c meson can only decay through weak 
interactions, it is an ideal system to study weak decays of heavy quarks. Due to the different 
decay rate of the two heavy quarks, the B c meson decays are rather different from those 
of B or B s meson. Although the phase space in c — > s transition is smaller than that in 
b — > c decays, but the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V CS \ ~ 1 is 
much larger than \V c b\ ~ 0.04. Thus the c-quark decays provide the dominant contribution 
(about 70%) to the decay width of B c meson At LHC experiments, around 5 x 10 10 B c 
events per year are expected 0, 0], which provide a very good platform to study various B c 
meson decay modes. 

In fact, the B c decays have been studied intensively by many authors j4-l|j||. In Ref.j3], for 
example, Dhir and Verma presented a detailed analysis of the B c decays in the Bauer-Stech- 
Wirbel (BSW) framework, while the authors of the Refs. studied the B c meson decays in 
the non-relativistic or relativistic quark model. In the perturbative QCD (pQCD) approach, 



furthermore, various B c decay modes have also been studied for example in Refs.|llHl3 



In this paper, we will study the semileptonic decays of B c — >• (nr., J/ ^f)Wi (here I stands 
for leptons e, ft, and r.) in the pQCD factorization approach [14j, |l5j. The lowest order 
diagrams for B c — > (rj c , J/fy)Wi are displayed in FigJH where B stands for B c meson and M 
for 7] c or J/\l/ meson, and the leptonic pairs come from the 6-quark's weak decay. In this 
work, we firstly calculate the g 2 -dependent form factors for B c — > (r/ c , J /^) transitions, and 
then we will give the branching ratios of the considered semileptonic decay modes. 

The structure of this paper is as below: after this introduction, we collect the distribution 
amplitudes of the B c meson and the r]c, J mesons in Sec. II. In Sec. Ill, based on the k^ 
factorization theorem, we calculate and present the expressions for the B c — > (r] c , J /^) tran- 
sition form factors in the large recoil regions. The numerical results and relevant discussions 
are given in Sec. IV, a short summary will also be included in this section. 



II. KINEMATICS AND THE WAVE FUNCTIONS 



In the B c meson rest frame, with the rriB c {m) stands for the mass of the B c (rj c or J/^) 
meson, and r = m/rriB c , the momenta of B c and rj c (J/ty) mesons are defined in the same 



way as in Ref. [16 



Pi = — ^(1, 1,0_l), p 2 = —^(rv^^rri ,0jj, 



(1) 



where t/ + = 7/ + a/t/ 2 — 1 and rj~ = rj — a/t/ 2 — 1 with the definition of the 77 as of the form 

m B . 



V 



2m 



1 + 



2 2 
m — q 

771d 



(2) 



where q 2 = [p% —P2) 2 is the invariant mass of the lepton pairs. The momenta of the spectator 
quarks in B c and rj c (J /^) mesons are parameterized as 



ki = (0,xi^=pfcij_) 



m Bc , m Bc _ 
k 2 = [x 2 -j=-rr] + , x 2 -j=-rv ,k 2± ). 



(3) 
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FIG. 1. The typical Feynman diagrams for the semileptonic decays B c — > (t] c , J/^f)lv, where B 
stands for B c meson and M for r\ c or J/^f meson. 



For the J/\l/ meson, we define its polarization vector e as 

e L =-L( rj +,-r 1 -,0 A _), e T = (0,0,1), 



(4) 



where and €t denotes the longitudinal and transverse polarization of the J/\I/ meson. 

In the calculations, one can ignore the kj< contributions of B c meson. Furthermore, one 
can assume that the b and c quark in B c meson are on the mass shell approximately and 
treat its wave function as a S function. In this work, we use the same distribution amplitude 
for B c meson as those used in Refs.[HI, H 



:[O + m Bc )7 5 0,B c (x)W 



with 



2V2N, 



:S(x - m c /m Bc ) 



where m c is the mass of c-quark. 

For pseudoscalar meson r) c , the wave function is the form of 



(5) 



(6) 



V2N ( 



=75 ^"(x) + m v J s (x)}. 



(7) 



The twist-2 and twist-3 asymptotic distribution amplitudes, <f> v (x) and <p s (x), can be written 
as 0,0 



f 

(j?(x) = 9.58 r J ^ x(l - x) 



[x 



1.97 



2y/W c 

fr)c 



x(l — x) 



1 0.7 



x(l — x) 



1 - 2.8x(l -x) 

0.7 



1 -2.8x(l -x) 



2y/2N c 

As for the vector J/\I/ meson, we take the wave function as follows, 

1 



® t j/M 



mj/^ L (j) L {x) +^ L 0(x) 



V2N, 



Tpq> T (x) 



(8) 



(9) 



And the asymptotic distribution amplitudes of J/^> meson read as 17 



Y K J Y K J 2^2N C 
</>Hx) = 10.94-^^(1 - 2x 



x(l — x) 



x(l — x) 



~l 0.7 



1 - 2.8x(l - x) 

~i 0.7 



x(l — x) 



<j) v (x) = 1.67-^=[1 + C2x - l) 2 



1 - 2.8x(l -x) 



1 -2.8x(l -x) 



0.7 



(10) 



Here, <p L and T denote for the twist-2 distribution amplitudes, and 0* and <j) v for the twist-3 
distribution amplitudes. 



III. FORM FACTORS AND SEMILEPTONIC DECAYS 



As is well-known, the form factors of Bt s \ — > P, V, transitions have been calculated 
by many authors in the pQCD factorization approach and other methods 19|, [20J], and the 
pQCD predictions are generally consistent with those from other methods. 

The B c — > r] c form factors induced by vector currents are defined in Refs. 



(vM\c(ohA°)\ B c(pi)) 



(Pl + ^) M - 



m, 



m 



W) + 



2^ , m B c 



m 



q^F (q 2 lll) 



where q = p\ — p 2 is the momentum transfer to the lepton pairs, and m is the mass of r\ c 
meson. In order to cancel the poles at q 2 = 0, F + (0) should be equal to F o (0). For the sake 
of the calculation, it is convenient to define the auxiliary form factors fi(q 2 ) and f2{q 2 ) 



(Vc(P2)\c(Qh»b(0)\MPi)) = /i(9 2 W + Mf/W 
In terms of fi(q 2 ) and f 2 (q 2 ) the form factor F + (q 2 ) and F (q 2 ) read 

1 



(12) 



F + (q 2 ) = ~ [fi(q 2 ) + f 2 (q' 



1 + 



q 



m 



13, 



1 - 



13, 



(13) 



The form factors F + and F (or f\ and f 2 ) of the B c — > r\ c transition are dominated by a 
single gluon exchange in the leading-order and in the large recoil regions. And in the hard- 
scattering kernel, the transverse momentum is retained to regulate the endpoint singularity. 
The factorization formula for the B — > rj c form factors in pQCD approach is written as Til 



(vMl c(0)^b(0)\B c ( Pl )) = g 2 C F N c / d Xl dx 2 d 2 k 1T d 2 k 2T 



dz + d 2 zi dy d 2 yy 
(2rr) 3 (2tt) 3 



xe^^(^(p 2 )|c;(y) C ,(0)|0)e^-(0|c Q (0)6' 5 (^|5 c ( P 0)^r d - (14) 



.'ik 



Here P, V, S denote the pseudoscalar, vector and scalar meson respectively. 
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In the transverse configuration b-space and by including the Sudakov form factors and 
the threshold resummation effects, we obtain the B — >■ 7] c form factors f\{q 2 ) and f 2 {q 2 ) as 
following, 



/i(<? 2 ) = 8Tnji 2 Bc rC F j dx\dx 2 j bidbib 2 db 2 4>b c ( x i) 

x {2 [(p s (x 2 ) - rx 2 (f) v (x 2 )} ■ hi(x 1 ,x 2 ,bi,b 2 ) ■ a s (ti)exp [S a b(ti)] 

xii] + (r] + (f) v (x 2 ) - 2<f) s (x 2 )) 



+ 



4r c (j) s (x 2 ) - 2r(f) v (x 2 ) + 



x h 2 (x u x 2 , 61, 6 2 ) • a s (t 2 ) exp [S ab (t 2 )] 



(15) 



f 2 (q 2 ) = &ixm 2 B C F j dx x dx 2 j b x db x b 2 db 2 4>b c ( x i) 



x |[2^(x 2 ) -4rx 2 (0 s (x 2 ) -#"(x 2 ))] 
x /ii(xi,x 2 ,&i,&2) • a s (ti)exp [-^(ti)] 

xi(r/V(^)-2^(x 2 )) 



+ 



4r0 s (x 2 ) -2r c v (a;2) 



"q" — 1 

x /i 2 (xi, x 2 , 61, 6 2 ) • asfe) exp [-S ab (t 2 )} 



(16) 



where CV = 4/3isa color factor, r is the same as in Eqs. (fT|3~]) . while r c = m c /m Bc , and m c 
is the mass of c-quark. The functions hi and h 2 , the scales t\, t 2 and the Sudakov factors 
S a b are given in the Appendix A of this paper. 

For the charged current B c — y nJ &i, the quark level transitions are the b — > cl&i transition 
with the effective Hamiltonian l23| 



(17) 



G 

U ef f(b -»■ clui) = -^V cb C7 M (1 - 75)6 • *7 M (1 - 7s)^, 
a/2 



where Gp = 1.16637 x 10 5 GeV 2 is the Fermi-coupling constant. With the two form 
factors F + (q 2 ) and F (q 2 ), we can write down the differential decay width of the decay mode 
B c -> TjM as @, H 



rfr(E c -» r/ c /^) 
c?g 2 



G^^l 2 g 2 -m 2 



2 ) 2 ./( 



ml 



m* 



1927T 3 



m. 



2\2 



x{K + 2g 2 ) [g 2 - (m Bc - mf] [q 2 
+3m 2 (m| c -m 2 ) 2 F 2 (g 2 )}, 



4g2 

(m Bc +m) 2 ]F 2 (g 2 ) 



(18) 



where mi and m is the mass of the lepton and r] c respectively. If the produced lepton is e ± 
or /r^ the corresponding mass terms of the lepton could be neglected. 
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The form factors V(q 2 ) and Ao^ y2 (q 2 ) for B c — > J/^ transition are defined in the same 
way as in Refs. [l6|, Hi], 22] and are written explicitly as, 



V{q 2 ) = 47rm| c CV J dx\dx 2 J b\db\b 2 db 2 4>b c ( x i) ' (1 + r ) 

x |2 [4> T {x 2 ) - rx 2 (p v (x 2 )] ■ h 1 (xi,x 2 ,b 1 ,b 2 ) ■ a s (ti)exp [Sabih)] 



+ 



2r + 



Xi 



Vv 2 



<p V {x 2 ) 



h 2 (xi,X2,b 1 ,b 2 ) ■ a s (t 2 )exp [-Sab(t 2 ))} , (19) 



A (q 2 ) = 8nm 2 3c CF J dx\dx 2 J b\db\b 2 db 2 4>b c {x\) 

x | [(1 - r 2 x 2 + 2rx 2 r]) <p L (x 2 ) + r (1 - 2x 2 ) <p\x 2 )] 
x h x {x u x 2 , b u b 2 ) ■ a s (ti) exp [—Sabih)] 



+ 



2 , , xi xijrj + r{l - 2t] 2 )) . , 
r +r c + — ~ rxir] H ] <p (x 2 ) 



2^f^l 

h 2 (x 1 ,x 2 ,b 1 ,b 2 ) • a s (t 2 )exp [-S^fo)]}, 



(20) 



Ai(q 2 ) = 8tc7ti 2 Bc Cf / dx\dx 2 I b\db^b 2 db 2 4>b c { x i) 



1 + r 



x | [2(1 + rx 2 r])4> v (x 2 ) - 2(2rx 2 - r])<p T (x 2 )\ 
x h 1 (x 1 ,x 2 ,b 1 ,b 2 ) ■ a t (ti) exp[-Sab(ti)] 

+ [(2r c -xi + 2rr]) 4> v (x 2 )] ■ h 2 (x 1 , x 2 , b x , b 2 ) ■ a s (t 2 ) exp[-S ab (t 2 )] J, 



(21) 



A2(g2)= (l+r) 2 (^-r) 



Ai(q 2 ) — 8irm 2 Bc CF J dx\dx 2 J b\db^b 2 db 2 4>b c { x \) 



1 + r 



x { [(77(1 - r 2 a; 2 ) - r(l + a; 2 - 2x 2 ?7 2 ))</> L (:r2) + (l + 2r 2 x 2 - 7-77(1 + 2x 2 )) <£*(z 2 )] 
x x 2 , 61, 6 2 ) • exp [~S ab (ti)} 

1 O' [..:,) 



+ 



xi \rr] - - 



. > Xl \ I •X'l •■> 

l+(r c + r - y J 7/ + r (1 - r c - y + 



x h 2 {xi,x 2 , b u b 2 ) ■ a s {t 2 ) exp [-S ah {t 2 )\ 



(22) 



wherte r = mj/^/m^, Cp and r c are the same as in Eqs. fll5|16p . The expressions of the 
hard function h\ and h 2 , hard scales t\ and t 2 , and Sudakov function are all given in 
the Appendix A. One should note that the pQCD predictions for the form factors f\. 2 (q 2 ), 
V(q 2 ) and A ,i, 2 (g 2 ) as given in Eqs. (jl5|16|19H22p are all leading order results. The NLO 
contributions to the form factors of B — >■ (n, K) transitions as given in Refs. 0,8 do not 
apply here because of the large mass of c-quark and (r) c , J/^f) mesons. 
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The effective Hamiltonian for the decay modes B c — > J j^lv\ is the same as B c — » r\ r lvi. 
but corresponding differential decay widths are different. For B c — > J/tyl&i, we havej9|,|27, 28 



dT L {B e ^J/mv l )_ G 2 F \V cb \ 2 q 2 - mf Uq 2 -m 2 f (m 



Br: 



m- 



-q 2 \ 



dq 2 



192vr 3 m| c (q 2 ) 2 



x { 3m{X(q 2 )A 2 (q 2 ) + 



Q 

m 2 + 2q 2 
4m 2 



m 2 - q 2 ){m Bc + m)A 1 (q 2 ) 



4q 2 



\{q 2 



-\ 2- 



m Br + m 



Aiiq 2 



; (23) 



dT ± {B c -> J/mvt) 
dq 2 



G 2 F \V ch \ 2 q 2 



where m = m j/y, and X(q 




combined transverse and total differential decay widths are defined as 



Am Bc m is the phase-space factor. The 



dT, 



dT, 



+ 



dT. 



dq 2 dq 2 dq 



dT 

dq 2 



dT r dTr 
— - H 

dq 2 dq 2 



(25) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



In the numerical calculations we use the following input parameters (here masses and 
decay constants in unit GeV ) [HI, |29[ 



A ifi 4) = 0-287 



m 



>lc 



2.981, 



mj/y = 3.097, m,B c = 6.277, 
m c = 1.275 ± 0.025, m T = 1.777, f Bc = 0.489, t Bc = (0.45 ± 0.04) ps, 
\V cb \ = (41.2±5:|) x 10- 3 , f Vc = (0.420 ± 0.050), f J/9 = (0.405 ± 0.014). (26) 

By using the expressions in Eqs. 



and the definitions in Eq. ( ITS]) we calculate 
the values of the form factors F (q 2 ), F + (q 2 ), V(q 2 ), A (q 2 ), Ai(q 2 ) and A 2 (q 2 ) for given 
value of q 2 in the region of < q 2 < {m Bc — m) 2 . But one should note that the pQCD 
predictions for the considered form factors are reliable only for small values of q 2 . For the 
form factors in the larger q 2 region, one has to make an extrapolation for them from the 
lower q 2 region to l arg er q 2 region. In this work we make the extrapolation by using the 
formula in Refs. |9L 31 



F(q 2 ) = F(0)-exp[a-q 2 + b-(q 2 ) 2 ]. 



(27) 



where F stands for the form factors F 0t+ , V, A ^ j2 , and a, b are the parameters to be deter- 
mined by the fitting procedure. 

The numerical values of the form factors F 0j+ ,V and ^0,1,2 at q 2 = and their fitted 
parameters are listed in Table HI The first error of the pQCD predictions for the form factors 
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FIG. 2. The pQCD predictions for the g 2 -dependence of the form factors Fq, F + , V, Aq, A\ and Ai- 
With the solid lines stand for the central values, and the bands show the errors of the corresponding 
form factors. 



in Table d comes from the uncertainty of the decay constants of r] c and/or J/\I/ mesons, and 
the second error comes from the uncertainty of m c = 1.275 ± 0.025 [29(]. For the parameters 
(a, b), their errors from the decay constants f Vc , fj/y or from m c = 1.275±0.025 are negligibly 
small and not shown here explicitly. As a comparison, we also present some results obtained 
by other authors based on different methods in Table [Til One should note that the definition 
of the B c — > J/ty transition form factors in this paper are different from those in Ref. 0] 
(ISK) and [E| (HNV). From Table HTl we find that our results agree well with the results in 
other literatures. In FigfS], we show the pQCD predictions for the g 2 -dependence of those 
form factors, where the solid lines stand for the central values, and the bands show the 
theoretical errors of the corresponding form factors. 



TABLE I. The pQCD predictions for form factors i*b,+j V, A), 1,2 at q 2 = and the parametrization 
constants a and b for B c — > r] c and B c J/^Sf transitions. 





F(0) a b 




0.48 ±0.06 ±0.01 0.037 0.0007 




0.48 ±0.06 ±0.01 0.055 0.0014 


yB c ^J/-f 


0.42 ±0.01 ±0.01 0.065 0.0015 


^0 


0.59 ±0.02 ±0.01 0.047 0.0017 


A B C ^J/^ 


0.46 ±0.02 ±0.01 0.038 0.0015 


A 2 


0.64 ±0.02 ±0.01 0.064 0.0041 



By using the relevant formulas and the input parameters as defined or given in previous 
sections, it is straightforward to calculate the branching ratios for all the considered decays. 



S 



TABLE II. B c — > rj c ,J/ty transition form factors at q 2 = evaluated in this paper and in other 
literatures. 





pQCD 


WSL[9] 


EFG[7] 


ISK[6] 


HNV[5] 


DV[4] 




0.48 


0.61 


0.47 


0.61 


0.49 


0.58 




0.42 


0.74 


0.49 


0.83 


0.61 


0.91 


A B c ->J/9 
^0 


0.59 


0.53 


0.40 


0.57 


0.45 


0.58 




0.46 


0.50 


0.50 


0.56 


0.49 


0.63 


,B c ^J/<t> 
A 2 


0.64 


0.44 


0.73 


0.54 


0.56 


0.74 



By making the numerical integration over the physical range of q 2 , we find the pQCD 
predictions for the branching ratios of considered decay modes: 

Br (B- Vc e-u e (fM-9,)) = (AAltUKfJ ± 0.39^)1^(^)1^^)) x 1(T 3 , 

bt(b- ^ Vc r-v T ) = (i.37jS;g(/*)±0.i2(7s o )jSS(y e6 )j2:SK)) * io~ 3 > 

Br {B- J/* e -j7 e (^-p M )) = (I0.03ig;g(/ J/V ) ± 0.89(r Bc )l°jl(y c6 )i^(m c )) x 10" 3 , 

fir (S- -> J/*r-v T ) = (2.92jJiS(/ J/w ) ±0.26(r Bc )l^ 6 7 (K fe )l^(m c )) x 10T 3 .(28) 

where the major theoretical errors come from the uncertainties of the input parameters /^ c , 
I Kb I) t b c and m c as given explicitly in Eq. (EBj) . 
From the pQCD predictions for the form factors -Fo,+, V, Ao,i,2 as given in Table I and the 
pQCD predictions for the branching ratios of B c — > (r] c , J/^)lv as given in Eq. fl28|) . we find 
the following points: 

(i) The form factor F (0) equals to -^+(0) by definition, but they have different q 2 - 
dependence. The error bands of Fo(q 2 ) and F + (q 2 ) in Figj2] are larger than that 
of V(q 2 ) and Ao,i t 2(q 2 )- The reason is that the uncertainty of the decay constant f Vc 
in B c — > 7] c transition is larger than the one of fj/q in B c — »• J/ty transition. 

(ii) The pQCD predictions for the form factors as listed in Table HT1 agree well with those 
obtained by using other methods. 

(iii) The pQCD predictions for the branching ratios of the four decay modes B c — > 
(j] c , J/ty)lv are at the order of 10~ 3 . Because of its large mass of r lepton, the 
decays involving a r in the final state has a smaller decay rate than those with light 
e~ or jjT . Since the ratio of the branching ratios has smaller theoretical error than the 
decay rates themselves, we here define two ratios R Vc and Rj/y, the pQCD predictions 
for them are the following 

ft- = f, ( ff - "''-"I " 3.2. for / = («,„), (29) 



Br (B~ -)■ r} c T-v T 
Br (B~ -> J/m~ 
Br (B~ JJWFu 

These relations will be tested by LHCb and the forthcoming Super-B experiments 



*>/• = f, ( ff - J JTl*!\ » 3.4, for l = ( e ,„). (30) 
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In short we calculated the branching ratios of the semileptonic decays B~ — > (r) c , J/i/;)l~&i 
in the pQCD factorization approach. We first calculated the relevant form factors by em- 
ploying the pQCD factorization approach, and then evaluated the branching ratios for all 
considered semileptonic B c decays. Based on the numerical results and the phenomenological 
analysis, we find that 

(i) For B c — > (i] c ,J/^/) transitions, the LO pQCD predictions for the form factors 
F 0j+ (0), ^(0) and ^0,1,2(0) agree with those derived by using other different meth- 
ods. 



(ii) The pQCD predictions for the branching ratios of the considered decay modes are: 



Br {B- -> Vc e-u e (fi-^)) = (4.4lj£gjj) x 10" 3 , 
Br (B- rj c T-Ur) = (1.37lg;ll) x lO" 3 , 
Br(B- -> J/^e-u e (fi-p„)) = (10.03+^ 3 ) x 10~ 3 , 

Br (B~ -> J/Vt-Vt) = (2.92±g:f2) x 10" 3 , (31) 



where the individual theoretical errors in Eq.( l28|) have been added in quadrature. 

(iii) We also defined two ratios of the branching ratios R Vc and Rj/m and presented the 
corresponding pQCD predictions, which will be tested by LHCb and the forthcoming 
Super-B experiments. 
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Appendix A: Relevant functions 

In this appendix, we present the functions appeared in the previous sections. The thresh- 
old resummation factors S t {x) is adopted from Ref. (l9j : 

and we here set the parameter c = 0.3. The hard functions h\ and h 2 come form the Fourier 
transform and can be written as 

h\{x\,xi,b\,b)) = K (Pibi)[9(bi - 62)^0 (ai6 2 )^o(ai 61) , 

+ 6(b 2 - hWaibJKoia^lStfa), [ ' 



h 2 (x 1 ,x 2 ,b 1 ,b 2 ) = #o(#2&2)[0(&i - b 2 )lQ{a 2 b 2 )KQ{a 2 b 1 ) 

+ 6(b 2 - 6i)/o(a2&i)^o(a 2 &2)]5 i (a; 1 ), { ' 



with ai\ = m Bc y/x 2 rri, f3\ = mB c y/xiX2rri + , a 2 = m Bc y/xirrj + and f3 2 — Pi- The functions 
K and I are modified Bessel functions. 
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The factor exp[—S a b(t)] contains the Sudakov logarithmic corrections and the renormal- 
ization group evolution effects of both the wave functions and the hard scattering amplitude 
with Sab(t) = S Bc (t) + S M (t), where 

S Bc (t) = s(x 1 ^,b 1 ) + 2 f ^ q (a s m, (A4) 

V 2 Ji/b 1 ^ 

S M (t) = s(x 2 ^,b 2 ) + s((l-x 2 )^,b 2 ) + 2 f ^ q (aM), (A5) 

V 2 V 2 Ji/b 2 ^ 

with the quark anomalous dimension 7 g = —a s /ir. The explicit expressions of the functions 
s(Q, b) can be found for example in Appendix A of Ref. [15] . The hard scales ti in Eqs.( TT6f22|) 
are chosen as the largest scale of the virtuality of the internal particles in the hard 6-quark 
decay diagram, 

t\ = max{m Bc y/x 2 rr], l/b 2 }, t 2 = max{m Bc v / a^rr/+, 1/bx, l/b 2 }. (A6) 
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